Symmetry of the superconducting order parameter in frustrated systems determined 

by the spatial anisotropy of spin correlations 
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We study the resonating valence bond (RVB) theory of the Hubbard-Heisenberg model on the 
half- filled anisotropic triangular lattice (ATL). Varying the frustration changes the wavevector of 
maximum spin correlation in the Mott insulating phase. This, in turn, changes the symmetry of the 
superconducting state, that occurs at the boundary of the Mott insulating phase. We propose that 
this physics is realised in several families of quasi-two-dimensional organic superconductors. 

PACS numbers: 



One of the major themes in condensed matter physics 
over the last few decades has been the deep connection 
between magnetism and unconventional superconductiv- 
ity. This is one of the key ideas that has emerged from the 
study of the cuprates [if, ruthenates 0] , cobaltates 0, 4 1 , 
heavy fermions [5|, organic superconductors @, 3 He 7|, 
and ferromagnetic superconductors Q. In the cuprates 
d x 2_ y 2 symmetry superconductivity emerges from the 
dopingof a Mott insulator with Neel order. Many the- 
ories [H, S, EH E3, includ 

ing RVB, suggest that in the 
metallic state spin correlations which are maximal near 
the wavevector (tt, tt) mediate superconductivity. In RVB 
theory [l|, Q superconductivity arises from the same 
strong correlations that give rise to antiferromagnetism in 
the Mott insulator. Alternative theories of the cuprates 
emphasize instead the role of different physics, such as 
stripes, phase fluctuations, or orbital currents 11 1. 

When frustration is introduced into a system the insu- 
lating state may not be Neel ordered and the spin cor- 
relations may not be strongest at (ir,ir). Therefore, a 
natural question to ask is what kinds of superconducting 
states do we expect to find when the spin correlations 
are different from the commensurate (jr, tt) correlations? 
In this Letter we study an RVB theory of the Hubbard- 
Heisenberg model on the half filled ATL to investigate 
this question. This is partially motivated by the fact 
that this model may describe whole families of organic 
superconductors . We find that as we vary the frustra- 
tion in our model the peak in the spin fluctuations in the 
insulating state moves continuously from (tt, tt), charac- 
teristic of the square lattice, via (2tt/3, 2tt/3), character- 
istic of the triangular lattice, to (n/2, 7r/2) characteristic 
of quasi-one-dimcnsional (qld) behaviour. This changes 
the symmetry of the superconducting state (see Fig. []}, 
from '<i x 2_j / 2' for weak frustration to 'd + id' at the maxi- 
mum frustration to l d xy + s' in the qld regime. We argue 
that these effects are realised in organic superconduc- 
tors such as K-(ET) 2 Cu[N(CN) 2 ]Br, k-(ET) 2 Cu 2 (CN) 3 , 
/3'-[Pd(dmit) 2 ] 2 X, and /3-(ET) 2 I 3 Q. 

The Hamiltonian of the Hubbard-Heisenberg model is 
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FIG. 1: (Color online.) Phase diagram of the Hubbard- 
Heisenberg model on the half filled ATL. As the frustration, 
J'/ J = (t'/t) 2 , is varied the spin correlations change from 
commensurate (n, n) characteristic of the square lattice for 
t'/t < 0.93, to incommensurate (q, q) in the highly frustrated 
regime, to commensurate (n/2,n/2) for t'/t > 1.3 character- 
istic of weakly coupled chains (see Fig. [2}. The spin correla- 
tions mediate superconductivity, and the changes in the spin 
correlations cause changes in the symmetry of the supercon- 
ducting state, which changes from l d x a_ y a' (A2) for small t'/t 
to 'd + id' (A 1 +iA 2 ) for t' ~ t to l s + d xy ' (Ai) for large t'/t. 



nihilates (creates) an electron on site i with spin a, Sj is 
the Heisenberg spin operator, and {ij} and (ij) indicate 
sums over nearest and next nearest neighbours across 
one diagonal respectively [H E3 (Fig. H]) . We study this 
model at exactly half-filling as this is appropriate for the 
/3, 0', k and A phase organic superconductors [6|. How- 
ever, studies of related doped models [§| suggest that the 
superconducting state evolves continuously upon doping. 

We study this Hamiltonian via the RVB variational 
ansatz |,E1E1, \RVB) =P G \BCS), where \BCS) is 
the BCS wavefunction and Pq is the partial Gutzwillcr 
projector, which we treat in the Gutzwiller approxima- 
tion. The problem reduces to solving the BCS and 
Gutzwiller variational problems simultaneously. This re- 
quires two mean-fields: Xk = 53 k , Vk— k'(cj^/tCk't) an d 



= J2k> Vk-k' (ck'TC-k'i)j where c k(7 is the Fourier 
transform of Ci CT , and d is the fraction of doubly oc- 
cupied sites. The pairing interaction, Vk = —6(1 — 
2e?) 2 [J(cos k x + cosfc y ) + J' cos(k x + k y )], arises from 
superexchange between nearest neighbours along the 
square (first term) an along one diagonal (second term), 
see Fig. [2j This potential directly links the symme- 
try of the superconductivity with the magnetic degrees 
of freedom. We assume singlet superconductivity. It 
then follows from the functional form of Vk and ba- 
sic trigonometry that the mean fields may be written 
as Ak = A x cos k x + A y cosk v + A ( jcos(fc a; + k y ) and 
Xk = Xx cos k x + x y cos k y + \d cos ( k x + k y ) - p, , the renor- 
malised chemical potential jl ensures half-filling. The 
symmetry of the ATL is represented by the group 
|6|. A basic theorem of quantum mechanics is that 
the eigenstates must transform like an irreducible rep- 
resentation of group which represents the symmetry of 
the Hamiltonian. It follows from this requirement that 
|A X | = |Aj,|, xl = X 2 y , and 8 = argA, = -argA y 0. 
Thus Ak = |A X | cos^cosfer +cos k y ) + \Ad\ cos(k x + k y ) 
+ i | A x | sin #(cos k x — cos k y ) . 

For simplicity we only consider Mott insulating states 
that are spin liquids, i.e., do not possess long-range mag- 
netic order. We are aware that for some parameters, e.g., 
large U/t and small t'/i, that states with magnetic or- 
der may have slightly lower energy. However, the d-wave 
spin liquid states considered here are quite competitive in 
energy [H,[23j]. Further, other work shows that the insta- 
bility of such ordered states to superconductivity, as U/t 
decreases, occurs for similar parameters as for spin liq- 
uid states. Hence, we suggest this simplifying assumption 
will not change our main results relating the supercon- 
ducting symmetry to the spatial anisotropy of the spin 
correlations in the parent Mott insulator. 

Within the Gutzwiller approximation, d = in the in- 
sulating phase and the model is equivalent to the Heisen- 
berg model. Therefore, results in the insulating phase do 
not explicitly depend on U . However, in the insulat- 
ing state the Hubbard model over-represents the Heisen- 
berg model. This leads to an SU(2) degeneracy of the 
insulating phase of the RVB theory [3,1§]. Physically 
this means that in the Mott insulator the mean fields 
are not physically distinct and the physical order pa- 
rameters are D — A 2 + x 2 = \J + X y and D' = 

yArf + Xd- It is straightforward to show that the spin 
correlations are peaked at the wavevector (q, q) where 
q = arccos (D 2 /2D' 2 ) . We solve the variational problem 
numerically on an 1000 x 1000 fc-space mesh. Fig. [^com- 
pares the wavevector found in this way from the RVB the- 
ory with the classical result, q — arccos (J/2 J') [l6j |. For 
J' I J < 0.87 (t'/t < 0.93) wc find that the spin correla- 
tions are commensurate and peaked at (tt,tt), consistent 
with a tendency towards Neel ordering. We also find 
commensurate spin correlations [peaked at (n/2,n/2)] 




FIG. 2: (Color online.) The wavevector, (q, q), where the spin 
correlations, which mediate superconductivity, are strongest. 
RVB theory shows that quantum effects increase the region 
where commensurate correlations are found relative to the 
classical theory. The lower inset shows the angle 9 which 
determines the symmetry of the superconducting order pa- 
rameter. 9 = — 7r/2 implies A2 ('d x 2_ y 2') superconductiv- 
ity; 9 — — n implies A\ ('s + d xy ') superconductivity; and 
— 7r/2 < 9 < —tv implies an l Ai + iA-z ('d + id') state. Note, 
in particular, that 9 = — 2iv/3 for t' — t, independent of 
U . The upper inset is a sketch of the ATL indicating the 
relevant hopping integrals (exchange parameters) to nearest 
neighbours (solid lines) and across one diagonal (dashed line). 



for large J' /J. This is the classical ordering wavevec- 
tor for uncoupled chains. It is difficult to determine ex- 
actly when the correlations becomes commensurate, as 
there is a smooth crossover (see Fig. [2]). However, it is 
clear that q - ir/2 for J'/ J > 1.7 (t'/t > 1.3). This 
shows that quantum effects enhance the stability of the 
region with commensurate spin correlations compared to 
the classical result. This effect is also found by other the- 
oretical methods [l7j]. In the region 0.87 < J'/ J < 1.7 
(0.93 < t' /t < 1.3) the insulating state is characterised 
by incommensurate spin correlations (except at the high 
symmetry point t' = t, sec below). In the metallic phase 
the two mean fields are physically distinct. Xk varies 
smoothly as the frustration, t'/t is varied (Fig. [31), but, 
three distinct superconducting phases are observed as is 
indicated by the behaviour of Ak (Figs. [21 El and 2]). 

For small t'/t we find that A^ = and 8 = —ir/2, thus 
Ak = A a; (cosfc a ; — cos fey). This is the prototypical form 
for a ^d x -2- y i' superconductor. Formally, Ak transforms 
according to the A2 representation of Ci v . This is consis- 
tent with the fact that for small t'/t the spin correlations 
in the insulating state are peaked at (tt,tt) (c.f., Fig. 

For large t'/t we find that A^A^A^ 7^ and 8 = 
—tt/2. Thus the order parameter takes the form Ak = 
I A;,; I (cos + cosk y ) + | A d I (cos k x cos k y — smk x smk y ). 
The first three terms are usually referred to as 'extended 
s (xs)' order parameters as they transform according to 
the trivial representation, but may have accidental nodes. 



2 



0.5 



1.5 




-1.5 





— r - 


















1 





0.5 



1 1.5 
t'/t 



FIG. 3: (Color online.) Variation of the RVB mean-fields as 
the frustration is varied for U = lOt. Xx,Xy> an d Xd all vary 
smoothly, but two phase transitions occur between supercon- 
ducting states. The first occurs when Ad becomes finite. But 
the second involves the angle 9 (inset to Fig. [2}. A x , A y , and 
Ad become small at large t'/t due because the bandwidth 
W oc t' in this regime and so we are moving away from the 
Mott transition as t'/t increases for t' ^> t (c.f., Fig. [TJ. 



xy 



However, the fourth term would be referred to as a l d 
state on the square lattice. On the square lattice these 
'xs' and l d xy " states belong to different irreducible rep- 
resentations of C± v . We therefore refer to this state as 
the 's + d X y state. However, we stress that Ak trans- 
forms solely as the A\ representation of C 2l ,; a direct 
consequence of the lower symmetry of the ATL. In this 
regime the spin correlations are (nearly) commensurate 
at (7r/2,7r/2), as expected for weakly coupled chains. 
Thus it is these qld correlations that cause the super- 
conductor to take 's + d X y symmetry. 

For t' ~ t; A x ,A y , A d ^ and -tt < 9 < -tt/2. Thus 
Ak has a non-trivial complex phase and breaks time re- 
versal symmetry (TRS): this might be detected by muon 
spin relaxation experiments 0, fl9l . 120] ] . The real part is 
the same as Ak for large t'/t and transforms according 
to the A\ representation. The imaginary part takes the 
same form as Ak for small t'/t and transforms according 
to the Ai representation. We therefore refer to this state 
either as the A\ + 1A2 or 'd + id' state. In this regime we 
have competition between spin correlations characteristic 
of the square lattice, which promote t d x 2_ y 2' supercon- 
ductivity, and those along the diagonal which favour a 
l s + d X y state. The compromise between these frustrated 
interactions is the A\ + 1A2 state with broken TRS. Fig. 
|4] details how Ak varies with the frustration for t' ~ t. 

Exactly at t' = t the lattice becomes hexagonal and 
has Cq v symmetry. In the insulating phase we find com- 
mensurate spin fluctuations peaked at (27r/3, 27r/3). In 
the superconducting state 'xs' terms transform like the 
A\ representation of Cq v . However, the l d x 2_ y 2 ' and l d xy ' 
terms transform as the E2 representation. E2 is a two- 
dimensional representation spanned by the l d x 2_ v 2' and 



'd xy ' terms. Thus, 'd+id' states that transform as the E2 
representation are expected on symmetry grounds on the 
hexagonal lattice for ap prop riate values of the Ginzburg- 
Landau coefficients 3, lj- This has already led to the 
prediction of broken TRS in K-(ET) 2 Cu 2 (CN)3_and f3'- 
[Pd(dmit) 2 ] 2 A on phenomenological grounds [19]. For 
t' ~ t ^ t' the lattice is slightly distorted away from C§ v 
symmetry and this leads to either the 'd x 2_ y 2 ' (for t' < t) 
or 's + d X y (for t' > t) component giving a greater con- 
tribution to Ak- This is clearly seen in Figs. [2] and [4] 
Experiments on CS2CUCI4 [2l| and other calculations [lTj 
suggest that RVB underestimates the size of the region 
where q ~ 27r/3. As these frustrated spin fluctuations 
drive A\ + iA 2 superconductivity this suggests that RVB 
theory may underestimate the stability of this phase and 
the size of the region of the phase diagram (Fig. [T|) where 
A\ + 1A2 superconductivity occurs. 

We note that our phase diagram (Fig. [1]) differs in 
the region around t' ~ t, from that recently proposed by 



others 22]. These differences arise because those works 
did not consider the possibility of insulating states with 
incommensurate spin correlations or l d + id' supercon- 
ducting states, and so found a spin liquid state (with 
commensurate spin correlations) for t' ~ t and large U/t. 

There is significant evidence that RVB physics is en- 
hanced on frustrated lattices 0, [23| and that it is relevant 
to layered organic superconductors @, [3, 22 1. Further, 
the RVB theory predicts a pseudogap in the metallic state 
above the superconducting critical temperature. Below 
about 50 K such a pseudogap is suggested by NMR re- 
laxation rate and Knight shift data [a, [25| . Additionally 
the insulating state of k-(ET) 2 Cu2(CN) 3 is a spin liquid 
0, 0] ■ Thus these materials provide an testing ground 
for the ideas presented here. The band structures suggest 
that these materials span the parameter range where the 
different superconducting order parameters occur Q . For 
example, t' < t in k-(ET) 2 Cu[N(CN) 2 ]C1 and k-(ET) 2 - 
Cu[N(CN) 2 ]Br which suggests that they have l d x 2_ y 2' 
(A 2 ) order parameters, t' ~ t in k-(ET) 2 Cu 2 (CN) 3 and 
/3'-[Pd(dmit) 2 ] 2 A and we propose that they have l d + id' 
(A\ + 1A2) order parameters, and t' > t in /3-(ET) 2 l3 
which suggests that it has an 's + d X y (Ai) order param- 
eter Q. These results are consistent with our current 
knowledge of the superconducting states of these materi- 
als, but much controversy remains over the experimental 
situation 0, [l9j]. It has also been argued that the su- 
perconducting state of the doped triangular lattice com- 
pound Naa;Co0 2 • yH 2 is an RVB state with l d + id' 
pairing [3, 0] . Ferromagnetic fluctuations are strong in 
doped triangular lattice systems [26( and so the possi- 
bility of triplet superconductivity needs to be considered 
carefully in both doped and half-filled systems 0- 

We have studied the RVB theory of the Hubbard- 
Heisenberg model on the ATL. Varying the frustration 
t'/t changes the spatial anisotropy of the spin correla- 
tions, which vary from being peaked on (tt, tt) for small 
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FIG. 4: (Color online.) The fc-dependence of superconducting order parameter, A(k) as the frustration (t 1 /t) is varied. Each 
plot covers the first Brillouin zone and the solid lines denote the non-interacting Fermi surfaces. For t'/t < 0.92 we have a 
'd x i _ y 2' superconductor where Ak transforms as the A2 representation of C2V and has symmetry required nodes along the lines 
k x — iky. This superconducting state is driven by the strong spin correlations at wavevector (n,iv). For t'/t > 1.06 we find an 
's + d xy ' order parameter which transforms like the representation of C2v This superconducting state is favoured by the 
strong spin correlations near (n/2,n/2). We find that this state has nodes, although they are not required by symmetry and 
therefore their location is dependent on t'/t and U /t. For t' = t the lattice maps onto the hexagonal lattice and the ground 
state is a 'd + id' state which transforms as the E2 representation of Ca v . Intermediate states such as those found at t'/t = 0.95 
and 1.04 still show the effects of the strongly frustrated triangular spin correlations and form l Ai + 4^2' states. The 'Ai + 1A2 
states found for 0.92 < t'/t < 1.06 all break time reversal symmetry and are fully gapped. The results shown are for U = lOt. 



t'/t, to incommensurate fluctuations for t' ~ t [except at 
t' = t where the (3tt/2, 3ir/2) fluctuations are commen- 
surate], to being commensurate at (tt/2,tt/2) for large 
t'/t. This drives changes in the symmetry of the super- 
conducting state. We propose that, as i d x 2_ y 2' results 
from proximity to a Neel ordered state, so 'd + id' super- 
conductivity arises from proximity to a spiral state and 
! s + id X y superconductivity is driven by qld spin fluctu- 
ations. The generality of the connection between (jr, 71") 
spin correlations and 'd x i_ y i' 0, [ij [nj EH suggests that 
our results are valid beyond the Hamiltonian studied and 
the approximations used in this Letter. This clearly begs 
the question: which superconducting states are driven by 
proximity to other magnetic orderings? 
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